An Unstructured Mesh Convergent Reaction-Diffusion Master Equation for Reversible Reactions

The convergent reaction-diffusion master equation (CRDME) was recently developed to provide a lattice particle-based stochastic reaction-diffusion model that is a convergent approximation in the lattice spacing to an underlying spatially-continuous particle dynamics model. The CRDME was designed to be identical to the popular lattice reaction-diffusion master equation (RDME) model for systems with only linear reactions, while overcoming the RDME's loss of bimolecular reaction effects as the lattice spacing is taken to zero. In our original work we developed the CRDME to handle bimolecular association reactions on Cartesian grids. In this work we develop several extensions to the CRDME to facilitate the modeling of cellular processes within realistic biological domains. Foremost, we extend the CRDME to handle reversible bimolecular reactions on unstructured grids. Here we develop a generalized CRDME through discretization of the spatially continuous volume reactivity model, extending the CRDME to encompass a larger variety of particle-particle interactions. Finally, we conclude by examining several numerical examples to demonstrate the convergence and accuracy of the CRDME in approximating the volume reactivity model.Comment: 35 pages, 9 figures. Accepted, J. Comp. Phys. (2018

A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems

Stochastic reaction-diffusion models have become an important tool in studying how both noise in the chemical reaction process and the spatial movement of molecules influences the behavior of biological systems. There are two primary spatially-continuous models that have been used in recent studies: the diffusion limited reaction model of Smoluchowski, and a second approach popularized by Doi. Both models treat molecules as points undergoing Brownian motion. The former represents chemical reactions between two reactants through the use of reactive boundary conditions, with two molecules reacting instantly upon reaching a fixed separation (called the reaction-radius). The Doi model uses reaction potentials, whereby two molecules react with a fixed probability per unit time, $\lambda$, when separated by less than the reaction radius. In this work we study the rigorous relationship between the two models. For the special case of a protein diffusing to a fixed DNA binding site, we prove that the solution to the Doi model converges to the solution of the Smoluchowski model as $\lambda \to \infty$, with a rigorous $O(\lambda^{-1/2 + \epsilon})$ error bound (for any fixed $\epsilon > 0$). We investigate by numerical simulation, for biologically relevant parameter values, the difference between the solutions and associated reaction time statistics of the two models. As the reaction-radius is decreased, for sufficiently large but fixed values of $\lambda$, these differences are found to increase like the inverse of the binding radius.Comment: 21 pages, 3 Figures, Fixed typo in titl

Modeling genetic circuit behavior in transiently transfected mammalian cells

Binning cells by plasmid copy number is a common practice for analyzing transient transfection data. In many kinetic models of transfected cells, protein production rates are assumed to be proportional to plasmid copy number. The validity of this assumption in transiently transfected mammalian cells is not clear; models based on this assumption appear unable to reproduce experimental flow cytometry data robustly. We hypothesize that protein saturation at high plasmid copy number is a reason previous models break down and validate our hypothesis by comparing experimental data and a stochastic chemical kinetics model. The model demonstrates that there are multiple distinct physical mechanisms that can cause saturation. On the basis of these observations, we develop a novel minimal bin-dependent ODE model that assumes different parameters for protein production in cells with low versus high numbers of plasmids. Compared to a traditional Hill-function-based model, the bin-dependent model requires only one additional parameter, but fits flow cytometry input-output data for individual modules up to twice as accurately. By composing together models of individually fit modules, we use the bin-dependent model to predict the behavior of six cascades and three feed-forward circuits. The bin-dependent models are shown to provide more accurate predictions on average than corresponding (composed) Hill-function-based models and predictions of comparable accuracy to EQuIP, while still providing a minimal ODE-based model that should be easy to integrate as a subcomponent within larger differential equation circuit models. Our analysis also demonstrates that accounting for batch effects is important in developing accurate composed models.Accepted manuscrip

Uniform asymptotic approximation of diffusion to a small target: Generalized reaction models

The diffusion of a reactant to a binding target plays a key role in many biological processes. The reaction radius at which the reactant and target may interact is often a small parameter relative to the diameter of the domain in which the reactant diffuses. We develop uniform in time asymptotic expansions in the reaction radius of the full solution to the corresponding diffusion equations for two separate reactant-target interaction mechanisms: the Doi or volume reactivity model and the Smoluchowski-Collins-Kimball partial-absorption surface reactivity model. In the former, the reactant and target react with a fixed probability per unit time when within a specified separation. In the latter, upon reaching a fixed separation, they probabilistically react or the reactant reflects away from the target. Expansions of the solution to each model are constructed by projecting out the contribution of the first eigenvalue and eigenfunction to the solution of the diffusion equation and then developing matched asymptotic expansions in Laplace-transform space. Our approach offers an equivalent, but alternative, method to the pseudopotential approach we previously employed [Isaacson and Newby, Phys. Rev. E 88, 012820 (2013)PLEEE81539-375510.1103/PhysRevE.88.012820] for the simpler Smoluchowski pure-absorption reaction mechanism. We find that the resulting asymptotic expansions of the diffusion equation solutions are identical with the exception of one parameter: the diffusion-limited reaction rates of the Doi and partial-absorption models. This demonstrates that for biological systems in which the reaction radius is a small parameter, properly calibrated Doi and partial-absorption models may be functionally equivalent

Reactive Boundary Conditions as Limits of Interaction Potentials for Brownian and Langevin Dynamics

A popular approach to modeling bimolecular reactions between diffusing molecules is through the use of reactive boundary conditions. One common model is the Smoluchowski partial absorption condition, which uses a Robin boundary condition in the separation coordinate between two possible reactants. This boundary condition can be interpreted as an idealization of a reactive interaction potential model, in which a potential barrier must be surmounted before reactions can occur. In this work we show how the reactive boundary condition arises as the limit of an interaction potential encoding a steep barrier within a shrinking region in the particle separation, where molecules react instantly upon reaching the peak of the barrier. The limiting boundary condition is derived by the method of matched asymptotic expansions, and shown to depend critically on the relative rate of increase of the barrier height as the width of the potential is decreased. Limiting boundary conditions for the same interaction potential in both the overdamped Fokker-Planck equation (Brownian Dynamics), and the Kramers equation (Langevin Dynamics) are investigated. It is shown that different scalings are required in the two models to recover reactive boundary conditions that are consistent in the high friction limit where the Kramers equation solution converges to the solution of the Fokker-Planck equation.Comment: 23 pages, 2 figure

The influence of molecular reach and diffusivity on the efficacy of membrane-confined reactions

Signaling by surface receptors often relies on tethered reactions whereby an enzyme bound to the cytoplasmic tail of a receptor catalyzes reactions on substrates within reach. The overall length and stiffness of the receptor tail, the enzyme, and the substrate determine a biophysical parameter termed the molecular reach of the reaction. This parameter determines the probability that the receptor-tethered enzyme will contact the substrate in the volume proximal to the membrane when separated by different distances within the membrane plane. In this work, we develop particle-based stochastic reaction-diffusion models to study the interplay between molecular reach and diffusion. We find that increasing the molecular reach can increase reaction efficacy for slowly diffusing receptors, whereas for rapidly diffusing receptors, increasing molecular reach reduces reaction efficacy. In contrast, if reactions are forced to take place within the two-dimensional plasma membrane instead of the three-dimensional volume proximal to it or if molecules diffuse in three dimensions, increasing molecular reach increases reaction efficacy for all diffusivities. We show results in the context of immune checkpoint receptors (PD-1 dephosphorylating CD28), a standard opposing kinase-phosphatase reaction, and a minimal two-particle model. The work highlights the importance of the three-dimensional nature of many two-dimensional membrane-confined interactions, illustrating a role for molecular reach in control-ling biochemical reactions.Published versio